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Using the exact results of the 0(3) nonlinear sigma model (NLSM) and a few quantitative numer- 
ical data for integer-spin antiferromagnetic (AF) chains, we systematically estimate all magnon ex- 
citation energies of JV-leg integer-spin AF ladders and tubes in the weak-interchain-coupling regime. 
Our method is based on a first-order perturbation theory for the strength of the interchain coupling. 
It can deal with any kind of interchain interactions, in principle. We confirm that results of the 
perturbation theory are in good agreement with those of a quantum Monte Carlo simulation and 
with our recent study based on a saddle-point approximation of the NLSM [Phys. Rev. B 72, 
104438 (2005)]. Our theory further supports the existence of a Haldane (gapped) phase even in a 
d-dimensional (d > 2) spatially anisotropic integer-spin AF model, if the exchange coupling in one 
direction is sufficiently strong compared with those in all the other directions. The strategy in this 
paper is applicable to other JV-leg systems consisting of gapped chains which low-energy physics is 
exactly or quantitatively known. 
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I. INTRODUCTION 

JV-leg quantum spin systems, which we study in 
this paper, are a natural extension of a purely one- 
dimensional (ID) spin chain. The study of iV-leg sys- 
tems has continued for more than a decade. 1-5 Of course, 
it is generally more difficult to quantitatively solve and 
understand JV-leg systems with a larger leg number JV, 
although in (1+1)D systems, there are various powerful 
theoretical strategies (for example, conformal field the- 
ory, bosonizations, integrability methods, exact diagonal- 
ization, density- matrix renormalization group (DMRG), 
Monte Carlo methods, etc.). Actually, many theoretical 
studies for iV-leg systems have focused on two-leg sys- 
tems, while three- or higher-leg systems have not been 
thoroughly investigated. However, some magnets with a 
larger JV indeed exist and have been investigated exper- 
imentally. 6,7 Recently one can fabricate a few spin tube 
materials 8-13 as well as standard spin ladder ones: the 
former (latter) has a periodic (open) boundary condi- 
tion along the interchain (rung) direction. These facts 
motivate us to study spin systems with an arbitrary 
leg number. In addition, JV-leg systems pose interest- 
ing questions in the context of statistical physics: how 
do systems with a finite JV approach the correspond- 
ing higher-dimensional ones (infinite- JV systems), how do 
low-energy properties of JV-leg systems depend on the 
value of JV (for instance, an even-odd character exists or 
not), is a critical value of the strength of the interchain 
interaction finite or zero in an infinite-JV system, etc. 
These questions require systematic and quantitative pre- 
dictions for JV-leg systems. (The questions are partially 
resolved. 2 ' 14 ~ 22 ) 

Now, thanks to the above-mentioned theoretical tools 
in (1+1)D systems, some simple ID spin systems (Heisen- 



berg chains, two-leg ladders, etc.) and field theories (sine- 
Gordon model, nonlinear sigma models, etc.) have been 
exactly or quantitatively solved. Such results are not only 
important for the deep and quantitative understanding 
of the solvable systems themselves, but also very useful 
as a starting point to study more complicated or dec- 
orated systems such as large-JV ones. As well known, 
a low-energy effective theory for integer-spin antiferro- 
magnetic (AF) Heisenberg chains is the relativistic 0(3) 
nonlinear sigma model (NLSM) which is exactly solv- 
able. Furthermore, the low-energy physical quantities of 
integer-spin AF chains have been numerically estimated. 
In this paper, applying these established results of AF 
spin chains, we formulate the Rayleigh-Schrddinger-type 
perturbation theory for the interchain (rung) coupling in 
JV-leg integer-spin AF systems. As a result, all magnon 
dispersions are determined as a function of the leg num- 
ber JV, the strength of the rung coupling and the spin 
magnitude S, in the weak-rung-coupling regime. We ver- 
ify the validity of our perturbation theory by compar- 
ing results with those of quantum Monte Carlo (QMC) 
method. 

Here, we note that quite recently we have studied 
the same JV-leg integer-spin AF systems with a weak 
rung coupling, based on a different approach, the NLSM 
plus a saddle-point approximation (SPA) method. ' 
Although the SPA is not reliable for quantitative pre- 
dictions, we believe that the results provide the system- 
atic understanding of JV-leg integer-spin systems. The 
perturbation theory in this paper will also be useful in 
verifying the validity of the NLSM plus SPA. 

The outline of this paper is as follows. In Sec. II, 
we construct a perturbation theory for standard JV-leg 
integer-spin AF ladders and tubes within the lowest or- 
der of the rung coupling. It leads to an analytical formula 
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for all the magnon dispersions of the ladder and tube 
systems, elucidating several characteristics of low-energy 
excitation structure in the N-\eg systems. The predicted 
bands are indeed consistent with those of QMC method, 
and support our previous theory based on the NLSM 
plus SPA. In the next two sections, we apply pertur- 
bation theory to infinite-leg systems (higher-dimensional 
spin systems) , and a system with a generalized rung cou- 
pling. In the final section, a brief summary of this paper 
is presented, and the potential of perturbation theory is 
discussed. Appendixes A and B, respectively, provide the 
solutions of eigenvalue problems of simple matrices, and 
the spin correlation functions calculated by the NLSM 
framework. They are used in Sec. II. In Appendix C, we 
shortly discuss higher-order perturbation terms and the 
difficulty in them. 



II. PERTURBATION THEORY 

This section is the main content in this paper. The 
construction of perturbation theory and its basic results 
are presented. 

A. Model 

Throughout the paper, we mainly consider the fol- 
lowing Hamiltonian of standard iV-leg integer-spin AF 
Heisenberg ladders and tubes: 



(la) 



N 



H\\ = jJ2Y,Sn,j-S n , j+ l, J>0, (lb) 

n=l j 



N' 



T~i± — J± J]] S n j ■ S n +ij, 

n=l j 



(1c) 



where S n j is integer-spin operator on site j in nth chain. 
The first term is the Hamiltonian of independent N 
integer-spin AF chains. The second is the rung coupling 
term, namely the exchange interaction between neigh- 
boring chains. Ladders take N' = N — 1, while tubes do 
N' = N and SW+ij = Si,j- A remarkable point in this 
model is that if the leg number N is odd and the rung 
coupling Jj_ is positive (i.e., AF) in a tube, geometrical 
frustration along the rung exists. 



B. Integer-spin antiferromagnetic chains and 
nonlinear sigma model 

As stated in the Introduction, we will treat the rung 
coupling Ti.± as the perturbation term for the chain part 
Tin. For such a treatment, the exact or quantitative 



knowledge of the unperturbed part, the single integer- 
spin-5 1 AF Heisenberg chain H = JJ2j$j ' "j+ii i s 
necessary and essential. The chain H n has never been 
solved exactly, but its low-energy effective theory, the 
0(3) NLSM, is integrable and investigated well. We use 
the well-known structure of the NLSM in order to an- 
alyze the model (1). In this subsection, we review the 
structure and the relationship between the NLSM and 
the chain H . 23 ~ 27 

The action of the 0(3) NLSM is given by 



*NLSM 



dxdt 



- {d t rfi) 2 - v (d x m) 2 
v 



, (2) 



where rh(x, t) is the three-component vector field (x and 
t denote spatial and time coordinates, respectively), and 
the constraint m 2 = 1 is imposed. The coupling constant 
g and the velocity v should be determined as a function of 
parameters in the lattice system 7io- Within the approx- 
imation used in the original mapping by Haldane, 23-25 
v — > 2SJa (ao: lattice constant) and g — > 2/S, but true 
values of g and v somewhat deviate from these ones due 
to irrelevant-term effects neglected in the Haldane map- 
ping (see the discussion below and Table I) . In the NLSM 
picture, the original spin operator is written as 



Sj w (-l) x/a °Sm(x)+u(x), 



u = — m x o t m, 
gv 



(3a) 
(3b) 



where x — jao, and the uniform part u is the "angular 
momentum" of the NLSM, a Noether current (symbol 
x here denotes outer product). For the (1 + 1)D 0(3) 
NLSM, the vacuum docs not break any symmetries in 
the theory, and the Hilbert space consists of the vacuum 
and massive triplet particles. Let us introduce some no- 
tations. 26 ' 27 We can represent the vacuum and a onc- 
particle state with an 0(3) index a and a momentum p 
as follows: 



\p, a). 



(4) 



The index a (= 1, 2, 3) corresponds to the ath component 
of the field m. In the spin chain, these two states cor- 
respond to the singlet, unique ground state (disordered 
spin liquid) and an excited state with a spin-1 magnon 
around the momentum (wave number) k = n/a n . The 
energy dispersion of the triplet particle is given by 



(5) 



where A is the mass gap of the particle. This gap A 
should be regarded as the lowest excitation gap (i.e., 
Haldane gap) in the chain H . The true velocity v 
can be determined from the spin-1 magnon dispersion 
around k = ir/ciQ. Here, let us adopt the relativistic 
normalization for the states in Eq. (4), (0|0) = 1 and 
(P,a\p',f3) = 6 al3 4iruj{p)6{p - p')/v. From Eq. (4), a 
state with n triplet particles is expressed as 



\{p,a} n ) = \pi,ai;p2,a2;...;p n ,a n ), 



(6) 
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where the total energy and momentum are E n — 
S" = i and p n = Z)"=iPj> respectively. Using the 

notations in Eqs. (4) and (6), and the relativistic normal- 
ization, we can represent the resolution of the identity as 



1 



|0)(0| 



E n\ ^ 



n 



vdpj 



x|pi,ai; . . . ;p„,a n )(pi,ai; . . .;p n ,a n \ 

OC 

- E^ 

3=0 



(7) 



where Pjjto is the projection operator onto all the states 
with j particles, and V = |0)(0| [V] = Vj). 

The form-factor method and symmetry arguments tell 
us several matrix elements of the NLSM. A part of them, 
which will be used in the next subsection, is summarized 
below. 



(0\m a (x,t)\0) 


= o, 


(8a) 


(0\m a (x,t)\p 7 f3) 




(8b) 


(0\m a (x,t)\p,(3;p' 7l } 


= o, 


(8c) 


(p,P\m a (x,t)\p' n ) 


= o, 


(8d) 


(0\u a (x,t)\0) 


= o, 


(8e) 


(0\u a (x,t)\p,f3) 


= o, 


(8f) 


(0|u a (z,t)|p,/3;i/,7> 


= can be finite, 


(8g) 


(p,(3\u a (x,t)\p' 7l ) 


= can be finite. 


(8h) 



Equations (8a) (8c), (8d)and (8f) are consequences of the 
fact that one-particle states are odd for the symmetry 
operation m — ► —to, while the vacuum |0), two-particle 
states, and u are even. Since the vacuum is the sin- 
glet for the angular momentum u, Eq. (8e) holds (it is 
consistent with the fact that the ground state has zero 
magnetization in the chain 7io)- These symmetry argu- 
ments are also useful in analyzing more complicated ele- 
ments among multiparticle states. For the detail forms of 
Eqs. (8g) and (8h), see, for example, Refs. 28 and 29. The 
renormalization factor Z in Eq. (8b) could be determined 
in the theory space of the NLSM, as a function of g, v, 
and the ultraviolet cut off A. However, since the NLSM is 
considered as the effective theory of the spin chain in the 
present case, Z should be fixed so that the NLSM repro- 
duces low-energy properties of the chain. A proper way 
of fixing Z is given by the comparison of long-distance 
behavior of the spin correlation function calculated by 
the NLSM approach and that obtained numerically. We 
explain such a comparison method in Appendix B. 

Provided that A, v, and Z are accurately evaluated by 
numerical methods, the NLSM can quantitatively cover 
low-energy, long-distance properties of the integer-spin 
chain H.q. Fortunately, values of A, v, and Z have been 
indeed known in some cases. We summarize them in 
Table I. ft is known that Z is equal to g, 28 ' 29 if the 
NLSM is approximated as a model of triply degenerate 
massive free bosons via the SPA. Table I shows that the 



TABLE I: Data (A,v,Z) in spin-1, 2, and 3 AF chains Ho, 
calculated by some numerical methods (Refs. 28-31). 



Spin 


A 


V 


Z 


1 


0.410J 


2.49Ja 


1.26 


2 


0.0892J 


4.65Ja 




3 


0.0100J 







true value of Z in the spin-1 case is smaller than the 
Haldanc- mapping plus SPA value g — 2/S. This devia- 
tion must originate from the interaction among magnons. 
True values of Z and v would be closer to the Haldanc- 
mapping plus SPA ones, 2/S and 2SJao, respectively, 
with increasing S because the NLSM is a semiclassical 
approach for quantum spin systems. 



C. Lowest-order perturbation theory 

Let us develop a Rayleigh-Schrodinger-type perturba- 
tion theory in the A^-leg system (1), using the NLSM 
framework explained in Sec. II B. Being interested in 
low-energy properties of the model (1), we concentrate 
on the ground state and one-magnon states in the full 
Hilbert space of the unperturbed part 7i||. Following the 
notation in Sec. II B, we can represent the ground state 
and a one-magnon state as follows: 



m = ni°>«> 

n=l 

N 

\p,a,l) = [] l°>r 



(9) 

x\p,a)i, (10) 



where the index n (I) in |0)„ (\p, a)i) means the chain 
number, and thus \p, a, 1} stands for a state with one 
magnon in Ith chain. Similarly, the NLSM fields rh and 
u in lih chain are expressed as to; and ui, respectively. 
Here, we emphasize that all one-magnon states with a 
fixed momentum p are degenerate in the decoupled sys- 
tem 7i || , namely, the 3A^-fold degeneracy is present in the 
unperturbed one-magnon space. 

One can easily find that the first-order ground-state 
energy correction, Eq^ = ($|7ijJ$), is exactly zero, be- 
cause Eqs. (8a) and (8f) hold. This result E^g — is 
also satisfied in the original lattice system framework. 

We next consider one-magnon states. From Eq. (7), 
the projection operator onto all one-magnon states in N 
decoupled chains is given by 



= EE 



n— 1 a 



vdp 
Airoj{p) 



\p,a,n)(p,a,n\. (11) 



The effective Hamiltonian of one-magnon states in the 
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first-order perturbation is 



N' 



n'f = v one n ± v 



dx 
a 



(12) 



J±(S 2 m n ■ m n+1 + u n ■ u n +i) 



where uin+i = rh\, Un+i = Ui, and we dropped oscil- 
lating terms [oc (—l) x ' a °]. Using the matrix elements in 
Eq. (8), we can transform Wf 1 as 



■^lst 



EE 

N' 



vdp s 



vdpt , w , / S 2 J±_ 



N' 



■sr-^ v-^ / vdp S 2 Zv 

2^2^ I — — n — J± 



Airuj(p) 2uj(p)a 



\p, a,n)(p,a,n + l\ + \p,a,n+ 1) (p, a, n\ 
I 



(13) 



where \p, a, N + 1) = \p, a, 1) and we used the formula 
Y~ f dxe l ( p ~ p ) x = S(p —p')- The product u n ■ u n +\ docs 
not contribute to the first-order perturbation. One sees 
that the perturbation H± connects one-magnon states 
in neighboring chains. This effective Hamiltonian pro- 
vides the magnon-band splitting from the 37V-fold de- 
generate dispersion u)(p). The integrand in Wf 1 can be 
reexpressed as the N x N tridiagonal matrix, 



Ol(p) 



/ wl(p) 
w±(p) 



UJ 



\ 



'■■ uj±(p) 
w±(p) 



,(14) 



where w_i_(p) = S 2 ZJ±v/ (2uj(p)ao), a column and a row 
each denote a chain-number index, and thus lo' = lo±_ 
(= 0) in tubes (ladders). Appendix A shows that the 
eigenvalues of Q± are 



A r = 2cj±(p) cosfc r , where (15) 
7r r = 1, . . . , N ladders, 

,N — 1 (mod N) tubes. 



N+l ■ 

r = 0, 



For the tube interpreted as wave number along 

the rung direction. From these results, we can conclude 
that the 3./V-fold degenerate magnon band uj(p) is split 
into N sets of new triply degenerate ones, 



'(p) = LO(p) + X r (p) 



y/A 2 +p^V 



,2^,2 



J± cosk r , (16) 



by the first-order correction. Of course, the remaining 
triple degeneracy is attributed to the spin-1 triplet (a = 
1, 2, 3). The gap of each band Lu^ st (p) is given by 



A r = c^ st (0) = A + 



S 2 Zv 
Aa 



Ji cos k r . 



(17) 



Here, it is found that the new dispersion w^ st (p) is not a 
relativistic form, though the contributing perturbation 
term ^ n m n ■ rh n +\ is invariant under the "Lorentz" 
transformation (field m a is a Lorentz scalar). This con- 
tradiction is due to the fact that our perturbation theory 
does not cquivalcntly treat the energy and the momen- 
tum. However, it is resolved by the interpretation that 
the dispersion w^ st (jp) is the first-order expansion form of 
the completely relativistic dispersion, 



/ OC2 7„ 

UJ r {p) = J(A 2 +p 2 V 2 ) + — J ± cos k r + 0(Jl). (18) 



a 



This is the unique relativistic band with the expanded 
form, Eq. (16). Equations (16)-(18) are the main results 
of perturbation theory. 

From the discussion in this subsection, one would easily 
see that perturbation theory can treat any kind of rung 
couplings, and be applied to ./V-leg systems on a lattice 
with an arbitrary geometric structure besides ladders and 
tubes. 



D. Features of band splitting 

For the spin-1 case, that is most interesting and real- 
istic, all values of A, v, and Z have been obtained as in 
Table I. They lead to 



3.8 



y/1 + 37p 2 a 2 



J±. 



(19) 



Therefore, the magnon band structure for TV-leg spin- 
1 systems can be predicted in the perturbation theory 
scheme. Gaps A r of spin-1 systems with several N are 
summarized in Fig. 1. 

The new dispersion (16) and Fig. 1 elucidate several 
features of the magnon-band splitting. 

(i) For ladder systems (N' = N — 1), the band split- 
tings in both AF and ferromagnetic (FM) rung-coupling 
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cases occur in the same way. Namely, the band struc- 
ture depends on the magnitude of the rung coupling, 
\J±\, but does not on the sign of J±. This is because 
the term J± J2 n ™ n ' ^n+i changes its sign via the uni- 
tary transformation m„ =Gven — > — m„ =Gven , and the same 
transformation does not affect P onc and the effective the- 
ory for the unperturbed part, N decoupled NLSMs. The 
difference between AF- and FM-rung sides would appear 
in higher-order corrections of the term Jj_ u n ■ u n +i, 
which is invariant under the above transformation. 

(ii) For the tubes (N' — N), in addition to the magnon 
triplet, the extra double degeneracy w* st (p) = uj^(p) is 
present except for r = and r = y modes [note that 
cosfej. = cos(— k r )]. This degeneracy is attributed to the 
symmetry of the ir rotation with respect to the diameter 
of the cross section of tubes, i.e., the parity symmetry 
along the rung direction (see Fig. 3 in Ref. 19). There- 
fore, it must remain even if higher-order corrections are 
taken into account. 

(iii) Like the ladder case, in even-leg tubes, sign change 
J± Y^n ■ m n+ i — > —J± J2 n mn ■ f&n+i is possible via 
Tn n — even — > —rn n=cvcn . Therefore, the band structure in 
an even-leg AF-rung tube and that in the corresponding 
FM-rung one are identical. On the other hand, the band 
structure in an odd-leg AF-rung tube differs from that 
in the corresponding FM-rung one. This is because any 
transformation connecting AF- and FM-rung systems are 
absent in odd-leg tubes. 

(iv) In odd-leg tubes, the lowest excitation is the triply 
[sixfold] degenerate band WQ St (p) [ w ±(V-i)/2(p)] f° r the 
FM-rung [AF-rung] case: this could refer to an even-odd 
character in the tube system (1). The rung-coupling- 
driven lowest-gap reduction in odd-leg AF-rung tubes is 
smaller than that in odd-leg FM-rung ones. The asym- 
metric splitting between AF- and FM-rung sides in odd- 
leg tubes, and the even-odd character is due to geomet- 
rical frustration along the rung direction in odd-leg AF- 
rung tubes. 



gap 



gap of N-spin (single rung) 
problem 




J region of 
J»|J±| the perturbation fheory 

FIG. 2: Expected lowest-gap profile in the iV-leg integer-spin 
system (1). The symbol A means the Haldane gap of the spin- 
s' AF chain Tio- In the two- leg spin-1 case, the gap profile 
have already been predicted in Refs. 20, 32 and 33. 



(v) The lowest magnon band in FM-rung tubes (N > 
3) Wg 8t (p) and that in even- leg AF-rung tubes w^* 2 (p) are 
independent of the leg number N. On the other hand (as 
expected), that in all the other systems falls down with 
increasing N. 

(vi) Perturbation theory shows that the lowest excita- 
tion gaps in all the iV-leg systems are a monotonically 
decreasing function of |Jj_|. However, as in Fig. 2, from 
the consideration of the strong-rung-coupling regime, one 
sees that when — J± ( J±) is much larger than J, the low- 
est gap is close to the Haldane gap of the spin-A^ x S AF 
chain and does not vanish (the gap of the N spin problem 
in one rung, which order is Jj_). Predicting such a behav- 
ior would be beyond the scope of the present approach 
based on perturbation theory in the rung coupling. 

(vii) As we already stated in Sec. II B, it is expected 
that as S is increased, Z and v, respectively, become close 
to g = 2/S and 2SJa Q . This result suggests the relation 



X r (p) 



S: large 



4S 2 J 



v/A 2 + AS 2 J 2 p 2 a 2 



J±cosk r . (20) 



Thus, we can infer that the band-splitting width mono- 
tonically increases with the growth of S. 

The characteristic properties (i)-(vii) of the magnon- 
band structure in the weak-rung-coupling regime is com- 
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pletely consistent with our recent predictions based on 
the NLSM plus SPA analysis (see Sec. Ill A and Figs. 
4-8 in Ref. 19). The property (v) cannot explicitly be de- 
rived from NLSM plus SPA, but the right panel in Fig. 
7 of Ref. 19 implies it. We thus conclude that qualitative 
features of magnon bands predicted by our previous work 
in Ref. 19 are supported by the present study. For more 
detailed physical interpretation of the results (i)-(vii), see 
Sec. Ill A of Ref. 19. The present study based on per- 
turbation theory furthermore provides quantitative pre- 
dictions on band structures in the weak-rung-coupling 
regime. In particular, the analytical expressions as in 
Eqs. (16)-(20) cannot be obtained in the previous NLSM 
plus SPA approach. 



E. Quantitative comparison with quantum Monte 
Carlo Results 

Here, we compare quantitatively the lowest-excitation 
gaps of N-leg systems evaluated by perturbation theory 
and those by QMC method. The QMC method cannot 
directly calculate the excitation gaps. However, it may 
be extracted from the proper correlation functions that 
can be calculated by the QMC. 21 > 31 > 33 

In Table II, we summarize slopes of the lowest- 
excitation gaps of TV-leg spin-1 AF systems (1) in 
the weak-rung-coupling limit, i.e., — [dA m i n /d\ J±\] j ± ^o, 
where A m ; n is the minimal gap in {A r }. Here, the slope 
data of QMC method (presented by Matsumoto) is de- 
fined by (A min |j i=0 - A min ||j x | =0 . O ij)/(0.01J), where 
A m i n |j i= o is of course equal to the Haldane gap of the 
spin-1 AF Heisenberg chain H , A(S = 1) = 0.4105J. 
In addition to Table II, in Ref. 33, Todo and co-workers 
minutely investigate the slope of the gap of the two-leg 
spin-1 AF ladder in the extremely weak- rung-coupling 
regime, |Jj_| < 0.01J, using QMC method (see Fig. 9 
in Ref. 33). They show that the slope is ~ 3.72 in the 
AF-rung side, and is ~ 3.75 in the FM-rung one. This is 
very close to the slope predicted by perturbation theory, 
3.8. From results of Ref. 33 and the data in Table II, we 
can conclude that perturbation theory is quantitatively 
valid in the sufficiently weak-rung-coupling regime, irre- 
spective of leg number N. 

In order to estimate the slope numerically, it is nec- 
essary to obtain the excitation gap in high precision for 
small rung couplings. In QMC method, the precise es- 
timate of the gap is not easy because it has to be de- 
termined indirectly from correlation functions. On the 
other hand, our perturbation theory directly gives the 
slope, and thus is expected to be more reliable. Actually 
our prediction relies on the numerical estimates of A, v, 
and Z. However, these are quantities defined on single 
chains, and can be determined accurately by numerical 
calculations. Moreover, we stress that QMC method does 
not work out well for frustrated (odd-leg, AF-rung) tubes 
because of so-called minus-sign problem. The present ap- 
proach has an advantage that it can be applied even to 



TABLE II: Slope of lowest gaps in N-\eg spin-1 AF sys- 
tems (1) near the decoupling point J± = 0, namely, 
— [dA m i n /d| J±\]j ± ^(). The data of QMC is evaluated by 
Munehisa Matsumoto. 

slope of lowest magnon gap 



system 


rung 


QMC 


perturbation 


2-leg ladder 


AF 


3.5 


3.8 




FM 




3.8 


3-leg ladder 


AF 


5.0 


5.4 




FM 




5.4 


3-leg tube 


AF 




3.8 




FM 




7.65 


4-leg ladder 


AF 


6.16 


6.15 




FM 




6.15 


4- leg tube 


AF 


7.57 


7.65 




FM 




7.65 



theory 



those frustrated systems without any problem. 



III. N -> oo LIMIT AND HIGHER DIMENSIONS 

The comparison with QMC results presented in 
Sec. II E implies that first-order perturbation results are 
efficient even for large N (i.e., higher-order corrections 
are negligible at least in a sufficiently weak-rung-coupling 
case, even if N is large). Assuming this is valid even in 
the N —> oo limit, we can present a few predictions for 
spatially anisotropic d-D (d > 2) spin systems, from the 
results in Sec. II. 

In the N — > oo limit, the distribution of k r is dense. 
Therefore, the lowest-excitation modes w lst (p) in all the 
infinite- leg (i.e., 2D) systems take cosfc r — > —1 or +1. 
Particularly, for the FM-rung tube [even-leg AF-rung 
tube] cases, as we already mentioned in Sec. II C, the 
lowest mode cjg st (p) [u^* 2 (p)] always takes cosfco = 1 
[coskN/2 = — 1] regardless of N(> 3): namely, the lowest 
mode in the 2D system already exists in a corresponding 
ID one with arbitrary N. From these properties in the 
N — > oo case, one finds that the slope of the gap reduc- 
tion — [c£A m i n /o?| Jj_|] j ± -^o is a finite value 2uj±(0)/J± in 
the infinite- leg system (1). It means that [A] the Haldane 
phase (ID gapped spin liquid) still survives in the 2D spa- 
tially anisotropic integer-spin AF system if the exchange 
coupling in the chain direction is extremely strong, and 
[B] there is a finite critical value of Jj_ (JD which lies 
between the Haldane phase and the Neel ordered one. 
This prediction is consistent with results of previous stud- 
ies. 17 - 34 - 37 

For the 2D spin-1 spatially anisotropic AF Heisen- 
berg model [i.e., an infinite-leg spin-1 AF system (1) 
with J± > 0], the critical value has been numeri- 
cally evaluated: for instance, the chain-mean-ficld plus 
exact-diagonalization analysis, 34 the cluster-expansion 
method, 37 and the QMC simulation 17 conclude > 
0.025J, ~ 0.056J, and ~ 0.04365J, respectively. On 
the other hand, perturbation theory provides the slope 
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A = 0.41J 


A + 7.65 J l/ 


^\A-7.65 J± 


Jl / 




J 



FIG. 3: Lowest excitation gap and critical point Jl for the 
N — > co, spin-1 case predicted by perturbation theory. The 
quantity A is the Haldane gap of the spin-1 AF chain Ho- 
of the gap 2u>±(Q)/ J± ~ 7.65 for the spin-1, N — > oo 
case. As we explain in Fig. 3, a naive extrapolation of 
the linear behavior gives a vanishing of the gap at a finite 
Jj_. This may be identified with the critical value, within 
the present framework. It leads to |J£| ~ 0.05 J. This is 
again consistent with the above numerical estimates. We 
emphasize that the critical values Jj_ for both AF- and 
FM-rung cases have the same magnitude in the first-order 
perturbation theory. This prediction also agrees with a 
recent QMC result. 22 

These results and those in Sec. II indicate that per- 
turbation theory is quantitatively reliable for any A^-leg 
integer-spin AF system, including the N — > oo case. 

As one readily expects, the above discussion of the 2D 
system can be generalized to d-D lattice cases. Namely, 
following the similar calculation of the first-order pertur- 
bation in Sec. II C, we can show that for a A^ d_1 -leg sys- 
tem on a d-D lattice, the lowest-magnon-gap slope also 
remains finite at the N — ► oo limit, like the case of lad- 
ders or tubes. For instance, the lowest gap of a weakly 
coupled integer-spin AF chains on d-D hypercubic lattice 
is evaluated as 

a— 2 

where is the interchain exchange in ath direction 
(the first direction is that of the chain). Moreover, in 

the case of Jj_ = J± , the slope of the gap is determined 
as — dAd/d\J±\ = ^ d ~£a — '■ in the spin-1, 3D case, the 
slope is ~ 15.3, twice 2wj_(0)/ J±. Thus, our first-order 
perturbation theory suggests that a gapped phase exists 
if all the interchain couplings are weak enough in any d-D 



spatially anisotropic integer-spin AF system. This pre- 
diction is consistent with the general belief that gapped 
phases are robust against small perturbations. 



IV. RUNG-COUPLING DECORATIONS 

As we already mentioned, perturbation theory, in prin- 
ciple, can deal with any sorts of interchain interactions. 
Therefore, in this section, as a generalization of the stan- 
dard coupling TL± , let us briefly investigate the following 
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FIG. 4: Decorated rung couplings. 



typical frustrated rung coupling: 

TL' = TL± + TL1+TL2, 

"Hi = K\ S nt j ■ S n +ij+i, 

TL2 = K^S^Sn,] ' S n +i,j—i. (22) 

n,j 

This coupling is illustrated in Fig. 4. For example, an 
infinite- leg system with K 1 ^ and K2 = (K\ = K2) 
may be called a spatially anisotropic model on triangular 
lattice (spatially anisotropic Ji-J2-like model). It is note- 
worthy that in contrast to the standard case of TL' = TL± , 
when Hi( or 2) is present alone, odd-leg FM-rung tubes ex- 
hibits geometrical frustration. 

As in the case of TL' = TL±, the first-order correction 
of the ground-state energy, (<E>|7ii i2 | < I ) ), are zero. Fur- 
thermore, the effective Hamiltonians for TL± t 2 in the one- 
magnon space are easily obtained in the similar way from 
Eq. (12) to Eq. (13). The results are 



n\% = v onc n h2 v onc = E E / (S^a~ a Kl ' 2 ) i ]p > a ' n ) e±ipa °^ + i| + H - c ] ■ ( 23 ) 

The exponential factors exp(±zpao) come from the fields rh n {x ± do) in TL\^.- In ladder cases, if the rung coupling TL' 
contains only TL\ or TL2, the factors could be replaced with unity by taking the field redefinition m n (x±ao) — > rh n (x) 
at the intermediate stage in the derivation of Eq. (23). The replacement is physically reasonable, because a ladder 
system with TL' — TL±_ is equivalent to one with TL' = 7ii( or 2) and K\t OT 2) = J±- However, one has to note that 
such a replacement procedure is not allowed when more than two kinds of rung couplings are simultaneously present. 
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Using Eq. (23) and results of eigenvalue problems in Appendix A, we obtain the following magnon bands of AT- leg 
systems with a generalized rung coupling H': 

/& I I S2ZVa " 1 x f l^-^- ipa °-^e^°|co S (fc r + g) tubes, 

w r W - V ^ +PV + y/ At+jPv 2 I |J± - A^e^ -K 2 e 4pa °|cosfc r ladders, 1 j 



where the angle 6> is defined by J ± -K ie - rpa « -K 2 e ipa ° = 
\J± - Kxe'^ - K 2 e ipa °\e ie . In the dispersion relation, 
the phase of e ±lpa ° may be negligible around p = 0. 
While, when p becomes too large (order of ir/ao), and 
thus the deviation between e ±lpa ° and unity becomes 
clear, the present NLSM picture for the AT- leg system 
is less reliable. 

In the case of J± = K\ + K 2 , the magnon-band de- 
generacy around p = is not lifted by the first-order 
correction of the rung coupling H'. This must be due 
to frustration among H±, Hi, and H 2 . From this pre- 
diction, we can expect that the gapped-phase region in 
the infinite-leg system with a frustrated rung coupling 
are much larger than that in the standard system with 
H! = H± (see Sec. III). This is consistent with a recent 
work using DMRG. 38 



V. SUMMARY 

In this paper, utilizing only the information about 
integer-spin AF chains Ho (exact results of the NLSM, 
and values of A, v, and Z), we have formulated a 
first-order perturbation expansion in the rung-coupling 
strength for the A-leg integer-spin system (1). All the 
magnon excitation energies have been quantitatively pre- 
dicted as in Eqs. (16)-(18) and Fig. 1. Several features 
of the magnon excitation structure have been elucidated 
[see comments (i)-(vii) in Sec. II D]. The results of per- 
turbation theory are supported by the QMC data and 
agree with our recent study based on NLSM and SPA in 
Ref. 19. We further applied the perturbation method to 
cases of N — > oo and systems with a generalized rung 
coupling in Sees. Ill and IV, respectively. Particularly, 
in Sec. Ill, it is predicted that a Haldane state (gapped 
state with ID nature) survives in spatially anisotropic d- 
D (d > 2) integer-spin AF systems, if the exchange cou- 
pling of one direction are sufficiently larger than all others 
[see Eq. (21)]. This agrees with several previous studies 
based on different methods. 17,34-37 We expect that if fit- 
ting parameters A, v, and Z are more accurately calcu- 
lated, perturbation theory leads to almost exact predic- 
tions in the weak-rung-coupling regime. 

The perturbation theory approach has some advan- 
tages over other methods, (i) It can lead to not only 
the lowest magnon band, but also all other magnon ones. 

(ii) One can predict several quantities with analytical ex- 
pressions if nonuniversal parameters (A, v, Z) are given. 

(iii) Perturbation theory is expected to be applicable to 
all the systems with any leg number N . Therefore, as 



one found in Sees. II and III, the dimensional crossover 
behavior (how a ID AT-leg system approaches the cor- 
responding d-D one when N becomes large) can be par- 
tially described, (iv) It is possible to investigate any kind 
of interchain interactions. 

We have thus demonstrated that the very simple per- 
turbation theory can give nontrivial results on weakly 
coupled integer-spin AF chains. Within the first order, 
perturbation theory is reduced to an eigenvalue problem 
in finite dimensions thanks to various selection rules [see 
Eq. (8)] . The precise numerical data known for the single 
chain then enables us to make quantitative predictions in 
the weak-rung-coupling regime. 

Similar approach should also be useful for more general 
quasi- ID systems consisting of weakly coupled gapped 
ID systems. When the ID system is gapless, pertur- 
bation theory is not quite easy. Usually, we have to rely 
on renormalization-group arguments, as in the case of 
the well-studied two-leg ladder consisting of spin-i AF 
chains. 
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APPENDIX A: TRIDIAGONAL MATRICES 

Here, we briefly summarize solutions of eigenvalue 
problems of the following N x N Hermitian matrices, 
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two-point spin correlation is represented as 



A = 



B 



( M A 2 
A\ Ay 



\A 2 
[By B 2 

B% By 



V 



A\\ 



■•■ A 2 

A* 2 A, J 



(Ala) 



\ 



'•• B 2 
Bl By J 



(Alb) 



The case in which A 2 and B 2 contain an imaginary part 
is not discussed well in literature. 

The eigenvalue problem for the matrix A (N > 3) is 
easily solved by the Fourier transformation method. The 
eigenvalues a r and the corresponding eigenvectors A r are 
given by 

2r 

a r = Ay + 2\A 2 \ cos(fc r + 6), k r = — n, (A2a) 



'(e^e 2 *',...^™*-), 



(A2b) 



where the angle 6 is defined by A 2 = \A 2 \e t9 , r — 
0, ...,N - 1 (mod A), and A 2 = 1. On the other 
hand, the matrix B can be diagonalized by assuming 
that the nth component of each eigenvector satisfies 
b n oc e~ m< ^ sin(fcn) where the angle (j> is the argument 
of B 2 (B 2 = \B 2 \e 1 ^) and k is a constant. As a result, 
the eigenvalues f3 r and the eigenvectors B r are 



(3 r = By + 2\B 2 \ cosk r , k r = 



N + l 



(A3a) 



N - 



- T (e~^ sin k r , sin 2k r , 



..., e - lAr<#, sinAfc r ), (A3b) 



where r = 1, . . . , N and Bl = 1. 



APPENDIX B: SPIN CORRELATION 
FUNCTIONS IN INTEGER-SPIN AF CHAINS 



lim (S?S%) « (-l)^S 2 (m a (x 3 )m a (x k )) 

\j-k\^oo 



+{u a (x j )u a (x k )), 



(Bl) 



where Xj(k) — j(k) x a . The term (u a m a ) disappears 
due to the symmetry of fh — > —m. Because m and u 
are one-magnon and two-magnon fields respectively, the 
most relevant part at \j — k\ — > oo is the first staggered 
term. Using the projection operators in Eq. (7), we esti- 
mate it as follows: 

(m a (x)m a (0)) = (m a {x)\j^V s 1 \m a {®) 

S 

= {m a (x)Vym a (0)) + ... 



Zv 



dp- 



2ir 



4?r J y/A 2 + p 2 v 2 
Z K (\x\/Z) + ... 



(B2) 



where ^ = v/A is the correlation length and Kq(x) is 
the modified Bessel function. If correct values of the 
spin-wave velocity v and the Haldane gap A are already 
known, the comparison between the spin correlation de- 
termined by Eqs. (Bl) and (B2), and the one done by a 
numerical method can fix the parameter Z. Actually, the 
value of Z in Table I is determined by such a method. 28 ' 29 
In the similar way from Eq. (Bl) to Eq. (B2), the long- 
(imaginary)time behavior of the equal-position two-point 
spin correlation function is evaluated as 



<S?(t)S?(0)> 



S 2 (m a (x,T)m a (x,0)} 

S 2 Zv f , e -^ T 
dp 



4-7T ./ r U)(p) 
At 



S 2 Z 



V8ttAt 



(B3) 



where r(= it) > is the imaginary time, Sf(r) 



e n ° T S?e- n ° T , m a (x,T) 



(x)e- HaT , and 



Hq (x) is the Hankel function of the first kind. The 
form (B3) is an expected result from Eq. (B2) (because 
the NLSM is a relativistic system). 



We demonstrate how the NLSM scheme derives the 
asymptotic form of spin correlation functions in the 
integer-spin AF chain Ti. - As shown below, comparing 
the resulting form and that estimated by a numerical 
method, the factor Z, etc. can be determined. 

As we said in the main text, the low-energy physics of 
the chain Tto is described by the NLSM (2). From the 
formula (3), the long-distance behavior of the equal-time 



APPENDIX C: HIGHER-ORDER TERMS 

We briefly consider higher than second-order correc- 
tions in the A- leg system (1). In principle, following the 
standard manner of perturbation theory, one would per- 
form calculations up to any order. However, as one will 
see below, there is a difficulty in extracting quantitative 
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predictions from higher-order terms within the perturba- 
tion theory framework of the main text. 

Generally, a second-order perturbation term is repre- 
sented as 

where \A) is an eigenstate of the unperturbed system 
7i|| with an energy Ea and a momentum Pa, and sim- 
ilarly \Q) stands for a Q-magnon state with an energy 
Eq = Y, q oj{p q ) and a momentum P Q = J2 q Pg (Pg is 
the momentum of each magnon). The matrix element 
(A\H±\Q) in the integrand would contain 

/' dr 
a 

where the factor I(A, Q) is independent of x. There- 
fore, noticing that J dxe^ p Q- p ^ x = 2ttS(Pq - P A ) and 
S(Pq — Pa) 2 ~ L x 5(Pq — Pa) where L is the chain 
length, we expect that a second-order term generally gen- 
erates a quantity of O(L). Actually, one can easily verify 
that in the second-order ground-state energy correction, 
the term originating from the lowest-energy intermediate 
states with two magnons is O(L). Under the assump- 
tion that one-magnon states still take the lowest excita- 
tion when a finite rung coupling is added in decoupled 



chains Tiu, all O(L) second-order perturbation terms for 
the ground state must completely cancel out those for 
one-magnon states at least in the thermodynamic limit 
(L — > oo). Only O(L ) terms, if they are present, can 
contribute to the second-order correction of the magnon- 
band splitting in the limit. 

These arguments appear to suggest that almost all the 
second-order perturbation terms are O(L). However, the 
factor I(A, Q) in Eq. (C2) can possess a finite-size cor- 
rection; for instance, the parameter Z = \ (0\m a (0)\p, a)\ 2 
might be expanded as Z = Z(L) = Z + Z 1 /L + Z 2 /L 2 + 
. . . (Z in Table I stands for Zq). Furthermore, v and A 
must also have a finite-size effect. Such finite-size cor- 
rections could generate an 0(L°) term. Unfortunately, 
as far as we know, finite-size correction terms have never 
been quantitatively estimated. Therefore, we cannot de- 
rive quantitative predictions from the second-order per- 
turbation expansion. 

The similar difficulty also emerges in calculations of 
higher than third-order corrections. Thus, we can con- 
clude that it is generally hard to quantitatively calculate 
higher-order terms in the scheme in this paper. Using the 
above expansion Z(L) = Zq + Z\/L + . . . , one might ob- 
tain some qualitative features of higher-order corrections 
in the magnon-band splitting. 
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